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Abstract. The discrepancy between the Pleiades cluster distance based on Hipparcos parallaxes and main sequence 
fitting is investigated on the basis of Stromgren photometry of F-type stars. Field stars with the same metallicity 
as the Pleiades have been selected from the mi index and a technique has been developed to locate the ZAMS 
of these field stars in color-magnitude diagrams based on the color/temperature indices b — y, v — y, and /3. 
Fitting the Pleiades to these ZAMS relations results in a distance modulus of 5.61±0.03 mag in contrast to the 
Hipparcos modulus of 5.36±0.06 mag. Hence, we cannot confirm the recent claim by Grenon (1999) that the 
distance problem is solved by adopting a low metallicity of the Pleiades ([Fe/H]piciados ~ —0.11) as determined 
from Geneva photometry. The metallicity sensitivity of the ZAMS determined by the field stars is investigated, 
and by combining this sensitivity in all three color/temperature indices b — y, v — y, and /3 we get a independent 
test of the Pleiades distance modulus which support our value of 5.61 mag. Furthermore, the field star sample 
used for the comparison is tested against theoretical isochrones of different ages to show that evolutionary effects 
in the field star sample are not biasing our distance modulus estimate significantly. Possible explanations of the 
Pleiades distance problem are discussed and it is suggested that the discrepancy in the derived moduli may be 
linked to a non-spherical shape of the cluster. 

Key words, open clusters and associations: individual: Pleiades - Hertzsprung-Russell (HR) and C-M diagrams - 
stars: distances - stars: evolution - stars: abundances 



1. Introduction 

The Hipparcos Space Astrometry Mission has provided 
accurate absolute trigonometric parallaxes for roughly 
120,000 stars, which are distributed all over the sky, and 
hence relatively accurate distance measurements for stars 
at a much larger distance than previous obtainable from 
ground based observations. This has given the opportu- 
nity to compare distances of several open clusters derived 
from direct trigonometric measurements with those de- 
rived from main sequence (MS) fitting. For some of the 
clusters there are discrepancies between the derived dis- 
tances, but in most cases the differences are within the 
estimated uncertainties. An exception is the Pleiades for 
which the distance modulus derived using the mean of the 
Hipparcos parallaxes is almost 0.3 mag smaller than that 
derived using the MS fitting technique. Possible explana- 
tions of this anomaly are: 

1. The errors of the MS fitting technique may be underes- 
timated. This could arise from difficulties in the tech- 



nique itself, or it could be due to serious errors in the 
adopted chemical composition of the Pleiades cluster. 
2. There may be systematic errors on small angular scales 
in the sky of the Hipparcos parallaxes which are un- 
derestimated. This could bias the inferred distance to 
clusters that only cover a small angular area in the sky. 

If these possibilities can be excluded we may have to draw 
the important conclusion that the theory of stellar struc- 
ture and evolution, is incomplete or in other words: The 
Vogt-Russell theorem that the location of a star in the 
Hertzsprung-Russell diagram is uniquely determined by 
its mass, age, and composition is violated. It is this as- 
pect of the Pleiades distance problem that makes it so 
interesting. 

The mean parallax of the Pleiades cluster inferred 
from the Hipparcos data ranges from 8.60±0.24 mas 
(Mer milliod et al. 1997) to 8.45±0.25 mas (van Leeuwen 
1999 ). These parallaxes correspond to a distance interval 



Send offprint requests to: D. Stello 



of 116±3 pc to 118±4 pc or a distance modulus interval 
of 5.33±0.06 mag to 5.37±0.07 mag. These distance mod- 
uli should be compared with those foun d from the MS 
fitting method. Pinsonneault et al. ( 1998 ) find a distance 
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modulus of 5.60±0.05 mag, based on an extensive multi 
color MS fitting analysis. They make use of several open 
clusters to check for different possible error sources, and 
both isochrones and an empirical Hyades MS are used as 
the zero point of the ZAMS. Pinsonneault et al. (1998) 
suggest that the discrepancy between the results from 
Hipparcos and the MS fitting method is due to spatial 
systematic errors on small angular scales in the Hipparcos 
data (Pinsonneault et al. 1998, Fig. 20) which are larger 



than expected (Lindegren |1988| , |1989| , |1997[ ) . From a com- 
parison of the Pleiades M S with those of the Hyades and 
a Persei clusters Eggen (1998) also concludes that the 
Hipparcos parallax distance of the Pleiades may be in er- 
ror by some 10 %. An investigation of the possible spatial 
systematic errors in the H ipparcos data is performed by 
Narayanan & Gould (|l999|) who confirm the distance mod- 
ulus of around 5.6 mag by an estimate of 5.58 mag, though 
with a fairly large error of ±0.18 mag. Their investigation 
is based on a variant of the moving cluster method to get 
distances for each individual cluster member. On the ba- 
sis of the stated spatia l syst ematic errors in the Hipparcos 
data Robichon et al. ( 199S| ) make an extensive investiga- 
tion of this problem, which makes use of the method of 
analyzing Hipparcos intermediate data described by van 
Leeuwen & Evans (1998). They recalculate the Hipparcos 
parallaxes and find a distance modulus of 5.36±0.06 mag 
for the Pleiades which is the value adopted in this paper 
for comparisons. 



Recently, Grenon (1999) has claimed that the Pleiades 
distance problem is solved by adopting a cluster metallic- 
ity of [Fe/H] Pleiades = —0.11 as determined from Geneva 
photometry instead of [Fe/H]pioiadcs — 0.0 as determined 
from high resolution spectroscopy (Boesgaard & Friel 
1990| ). In particular, Geneva photometry points to a large 



difference in the metallicity of Praesepe and Pleiades 
(A[Fe/H] = 0.28). To a large extent this explains the off- 
set between the two clust ers in various color-magnitude 
diagrams (Robichon et al. |2000|) . 

The purpose of the current study is to investigate the 
discrepancy between the Pleiades distance inferred from 
the Hipparcos mean parallax and from the MS fitting 
method by comparing Stromgren uvhy — (3 photometry 
of Pleiades F-type stars with field stars having the same 
metallicity as the Pleiades. We start out with a presen- 
tation of the Pleiades cluster and field star data followed 
by the calibration formulas and selection of stars used for 
the rest of the reductions. Then the MS fitting analysis is 
described including the fitting technique, which has been 
developed to locate the ZAMS of the field star sample, and 
how we test the age range of the field stars. Furthermore, a 
metallicity sensitivity analysis is presented and finally we 
discuss possible answers to the Pleiades problem including 
the suggestion that the discrepancy in the derived distance 
moduli could be a real effect caused by the structure of 
the cluster. 



2. MS fitting with Stromgren photometry 

The Pleiades, shifted to the distance inferred from the 
Hipparcos mean parallax, are compared to nearby field 
stars of the same metallicity in a color-magnitude dia- 
gram. This is done to see if the ZAMS defined by the field 
stars (hereafter denoted ZAMSpidd) does coincide with 
the Pleiades MS. Further, the distance modulus of the 
Pleiades is inferred from a best fit by the comparison of 
the Pleiades MS and the ZAMSpieid- This investigation 
also includes a sensitivity analysis of the ZAMSpicid locus 
to changes in [Fe/H]. The analysis is performed for F-type 
stars, and in three color/temperature indices b — y, v — y, 
and (3, to check for robustness and possible hidden errors 
in the MS fitting technique. 

The advantage of using field stars to define the zero 
point of the distance modulus is that the investigation 
will be purely empirical, and not affected by some possi- 
ble lacking ingredients in the theoretical calculations of the 
isochrones. So this investigation only relies on the Vogt- 
Russell theorem (Sect. The shape of the Pleiades MS 
turns out to be very similar to the ZAMSpioid shape and 
hence we avoid the problem of fitting to isochrones which 
do not match the cluster MS at every temperature inter- 
val. Using F-type stars gives the opportunity to correct 
for interstellar reddening and to estimate [Fe/H], based 
on the Stromgren indices and available empirical calibra- 
tion formulas. The assumed shape of the ZAMSpieid used 
for the comparison between the Pleiades MS and the field 
stars is found by a second-order robust least squares fit to 
the Pleiades stars. 

2.1. Data 

The uvby — (3 data for the Pleiades is taken from Table 
II of Grawford & Perry (|l97|) (hereafter CP76), which 
contains members classified as F-type stars by CP76. The 
mean error of one observation, was determined from the 
internal scatter in the measurements of all the Pleiades 
stars (Table I CP76), and were given as follows, cr(/3) — 
0.011, a{b ~ y) = 0.009, cr(mi) = 0.011, and a{ci) = 
0.012 (CP76). Taking into account that each star has been 
observed about 5 times we obtain the following typical 
mean errors of the mean indices of one star: cr(/3) = 0.005, 
cr(6 ~y) = 0.004, cr(mi) = 0.005, and cr(ci) = 0.005 mag. 

The magnitudes of the stars in CP76 were adopted 
from Johnson & MitcheU ( |1958D (based on the UBV 
system), so to get the Stromgren values {y transformed 
to V) the star magnitudes are taken from the WEBDA 
database (Mermilliod 200C). For about half of the stars 
the Stromgren value is not available, and the magnitudes 
are taken to be the average of the V values from UBV 
observations (also given in the WEBDA database). It is 
noted that for stars with both uvby and UBV photometry, 
the V magnitudes agree within ±0.02 mag. 

The field star data has been taken from a catalo gue of 
~ 30.000 stars observed in uvby - (i (Olsen |l999|) . This 
catalogue has been made by merging five published cata- 
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logues; all by E.H.Olsen. The sample used in this investi- 
gation (F-type stars) is based on three of these catalogues 
(Olsen |1983| , |198S| , |1994D . The overall RMS internal er- 
ror of one observation is (t{V) — 0.005, a{b — y) = 0.004, 
cr(mi) = 0.006, cr(ci) = 0.007, and a{l3) = 0.007. These 
errors are the conservative ones, in some of the catalogues 
they were in fact ~ 0.002 mag smaller, but the conserva- 
tive ones are adopted in this investigation. The majority of 
the stars were observed only once and a few 2 or 3 times, 
so the errors in the mean photometric values per star are 
not significantly different from the errors stated above. 

The absolute magnitudes are derived using the redden- 
ing corrected apparent magnitudes Vq ( transf ormed from 
y), and the Hipparcos parallaxes (ESA 1997 ). The men- 
tioned possible spatial systematic error in the Hipparcos 
parallaxes does not affect the locus of the ZAMSpidd, be- 
cause it can only have effects on small angular scales, and 
the field stars are distributed "randomly" all over the sky. 
Thus the error contributions from the parallaxes in the 
zero point for the ZAMSpieW is the global parallax error, 
which is less than 0.1 mas (Arenou et al. 1997j ). 

In addition to the internal errors for the Pleiades pho- 
tometry by CP76 and for the field stars by Olsen (1999) 
there could be systematic differences between the two 
sources. Especially the f3 index is critical, because the red- 
dening and hence the corrected color indices (5 — y)o and 
(^^ — y)o is determined from f3 (Sect. ^.2.1 ). It is, however, 
very unlikely that systematic errors in f3 could be large 
enough to explain the offset between the Pleiades and the 
field stars, which is of the order of 0.02 mag in /3 if we 
adopt the Hipparcos distance modulus of the Pleiades. In 
this connection we note that photometric observations of 
the P index is quite straightforward; no extinction cor- 
rection is needed and the transformation to the stan- 
dard system is linear without color terms. As discussed 
by Olsen ( 1983| ), systematic differences in f3 obtained with 
different telescopes and filter sets are 0.005 mag at most. 
Furthermore, we note that two of the Pleiades F-type 
stars f rom C P76 (Hz II 739 and 948) happen to be in 
Olsen ( |1999| ). The differences (01sen-CP76) are 0.007 and 
0.008, respectively. This does not point to any large sys- 
tematic errors, and a correction for this difference would 
in fact increase the offset of the Pleiades with respect to 
the field stars, but of course we cannot draw any strong 
conclusions from two stars only. 

2.2. Calibration and selection 
2.2.1. Reddening 

To correct for reddening, the color excess is determined 
for individual stars as E{b — y) = {b — y) — [b — y)o, where 
{b — y)o is found by an iterative calculation based on the 
empirical calibration given by Crawford (1975a): 



{b-y)Q = 0.222-)-l.llA/?-h2.7A/32- 
0.05&0 - (0.1 -t-3.6A/3)(5mo 



(1) 



where 

Sco = Co - co^zAMs(/3) , Smo = TOo3yadcs(/3) - mo , (2) 

Co = ci - 0.20 • E{b -y) , mo = mi + 0.30 • E{b - y) , (3) 

and A/3 ~ 2.72 — (3. The standard relation between /9, 
mo,Hyadcs, and cq^zams is found by interpolat ion be tween 



the data points given in Table I by Crawford (1975a). The 
iterative procedure uses the four Stromgren indices for 
each star, and an initial guess for the color excess as input 
data. The output (individual color excesses) was obtained 
when E{b — y) converged at the 0.0001 level. The expected 
error for E[b — y) is found as 



a[E{b-y)]^[a\b-y) + a'[{b-y)o\] 



1/2 



(4) 



where 



cr2[(5-y)o] = cr2(A/3)(l.ll + 5.4A/3 + 3.6<5mo)2+ (5) 
cr2((5mo)(0.1 + 3.6A/3)2 + Q.Qb^a^{5co), 

with cr(A/3) = cf{(3), a{SmQ) ~ o'(mi), and (t{Sco) ~ 
(t(ci). An estimate of this error is obtained by inserting 
the mean values of A/3 and Stuq from the samples. For the 
Pleiades stars it is a[E{b — y)]pioiadcs — 0.008 mag, while 
it for the field stars is cr[E{b — y)]picid — 0.010 mag. 



2.2.2. [Fe/H] calibration 

To be able to select nearby field stars with the same metal- 
licity as the Pleiades, [Fe/H] is calculated for the Pleiades 
using the empirical calibrations of Nissen (1981). The cal- 
ibration formula is as follows: 



[Fe/H] = -[10.5 + 50(13 - 2.626)](5mo + 0.12 



(6) 



where 6mo is like in Eq. (|^), and the constant of 0.12 is 
the adopted [Fe/H]Hyadcs- The adopted [Fe/H]Hyados does 
not affect the result of the comparison of the Pleiades with 
the field stars, because it is only relative. But the estimate 
of the absolute [Fe/H]pieiados value is of course depen- 
dent on the assumed [Fe/H]Hyades- The expected error in 
this calibration is obtained by a procedure similar to that 
presented in Sect. 2.2.1 , but this time (7{dmo) = cr(mo) 
is used. The results are tT([Fe/H])pieiades = 0.07 and 
(T([Fe/H])pieid = 0.07. 

2.2.3. Reduction step by step 

The criterion for being an F-type star is set to be 2.59 < 
(3 < 2.72, which is the /3 range used by Crawford ( 1975a| ) 
and Nissen ( |1981| ) for their calibration formulas (redden- 
ing and metallicity). 

To keep as many Pleiades members as possible, no 
stars from Table II (CP76) are rejected as a start ex- 
cept Hz II 948 which appear to be a non-member in both 
CP76 and the WEBDA database. For each member star 
the redd ening is calculated by the procedure described 
in Sect. 2.2.1, and their photometry measurements are 
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individually corrected. The mean reddening obtained is 
{E{b - 2/))pioiades = 0.031 ± 0.004 magQ, and the star- 
to-star RMS scatter is S[E{b — 2/)]piciadcs — 0.022 mag. 
Compared with the expected error of the E{h — y) de- 
termination {a[E{b — y)]pioiadcs = 0.008 mag), this indi- 
cates significant star-to-star reddening differences across 
the cluster. The mean color excess for the Pleiades ob- 
tained here is in quite good agreement with former results 
obtained from other investigations, which in general are 
in the range of about 0.03-0.04 mag (e.g. Pinsonneault et 
al. |1998| used E{b-y) = 0.7 x 0.04 mag ~ 0.03 mag; here 
using the relatio n between E{b — y) and E{B — V) from 
Crawford [l975b| ). 

The reddening corrected mo values together with the 
(3 observations are then used as input in Eq. (||), to 
get the Pleiades metallicity. The mean value derived is 
([Fe/H]) Pleiades — 0.01 ± 0.02^. This value is in very good 
agreement with spectroscopic results, which mostly come 
out with a near so lar m etallicity for the Pleiades (e.g. 
Boesgaard & Friel |1990D . The star-to-star RMS scatter 
is 5([Fe/H])piciados = 0.13, which is somewhat larger than 
the expected error (o-([Fe/H])pioiadcs — 0.07). 

From the catalogue of field stars, used in this investi- 
gation, there are 12658 stars which have 2.59 < (3 < 2.72 
(thus F-type stars), but 1194 stars of this group do not 
have Hipparcos parallaxes, so the absolute magnitude 
could not be derived, and they are therefore rejected. For 
each star in the remaining sample the reddening is cal- 
culated, as described in Sect. 2.2.1, and the photometry 



of every star is individually corrected. Due to statistical 
fluctuations and a low mean reddening of the field star 
sample {{E{b — y)) = 0.009) some stars turn out to have 
slightly negative values of E{b — y). In order to avoid any 
bias these negative values were not changed. Finally, the 
reddening corrected toq values and the (3 observations can 
be used as input data in Eq. , to derive the metallicity 
for every star. 

We choose a metallicity range of —0.10 < [Fe/H] < 
0.12 which is comparable to the Pleiades mean metallic- 
ity plus/minus a representative estimate of the metallicity 
scatter, and the number of F-type field stars remaining in 
this interval is 3389. The mean [Fe/H] of the remaining 
sample is not equal to ([Fe/H])piciadcs because the metal- 
licity profile of the original field star sample peaks around 
[Fe/H] = —0.15, thus [Fe/H] for the remaining sample is 
slightly shifted (by 0.01 dex) to a lower [Fe/H]. But since 
this is only half the error of the mean of the Pleiades 
[Fe/H] the effect is ignored. 

An additional selection of the field stars is made on the 
basis of the relative parallax error. If the My^ vs. (6 — y)o 
diagram is considered, the error in the absolute magnitude 
of the field stars is affected by the errors in Vqi — 
t/)o, and the distance modulus (through tt, the parallax). 
From the errors given in Sect. 2.1 and p. 2.1 , the error 
in Vo can be estimated as: (j{Vo) [cr^(V^) + (4.3cr[i?(6- 



y)]Fm\dfV''^ = 0.044 mag. The size of the effect on cj{Mva) 
from (j[{b — y)o] depends on the slope of the ZAMS in 
the color region of interest. A test plot is made to find 
the approximate slope, and it is found to be ^ 12. The 
effect from cr(7r)/7r, on the distance modulus, is found by 
differentiating the relation (m — M) — 51og(i) — 5, with 
respect to 7r~^, where (ra — M) is the distance modulus. 
This now leads to the following expression: 

{Mv,) ^a\Va) + {12 a[{b~y)o\f + [2.17 a{TT)/^f .(7) 

The optimized choice of the upper limit of (j{'k)/'k^ is 
when 2.17cr(7r)/7r - max{cr(Vb), 12 cr[(6 - y)o]} 0.12 
mag (where (j[{b — ?/)o]Ficid = 0.010 mag; Eq. (^)), which 
suggests a{n)/TT ^ 0.05. If a similar consideration is made 
in the Mvq vs. (y — j/)o or Mv^ vs. {3 diagram, the sug- 
gestion would be o'(7r)/7r ^ 0.04. To avoid different sam- 
ples of field stars in the three investigation parts, parallax 
measurements to a 5% accuracy are chosen, which left a 
sample of 782 stars. This seems to be the most reason- 
able choice, since the sample is large enough to make a 
clear definition of the ZAMS, and a lower error always is 
desirable. 



2.3. MS fitting analysis 

The MS fitting analysis is carried out in three diagrams 
Mvq vs. (6 — y)o, Mvg vs. (w — y)o, and Mvq vs. f3. Often 
only the (b — y)o case will be illustrated in the figures, but 
the results of the other diagrams will be given. Though 
(3 is not a color, all three diagrams will in the follow- 
ing be denoted color-magnitude diagrams, and the spec- 
ification "F-type" stars will be omitted, thus "the stars" 
or "all the stars" simply refers to the sample of F-type 
stars used in this investigation. A test plot of the Pleiades 
showed that there was a outlying star which was signifi- 
cantly cooler than the rest of the sample, and ~ 1 mag 
above the Pleiades MS (it has spectral type GO; Mendoza 
1956). The star is rejected as a likely non- member. The 
Pleiades mean color excess and metallicity did not change 
significantly if the GO star was included in the sample or 
not. The changes were only 0.002 mag in the average color 
excess, and ([Fe/H]) changed by 0.003 dex. 



With the errors stated in Sect. 2.1 and Sect. 2.2.1 the 



^ This value is based on all stars from Table II of CP76 except 
the likely non-members Hz II 739 and Hz II 948 (see Sect. O). 



estimated error per one star in Mv,-, for the field stars 
[Eq. (^], is in the range of 0.14-0.17 mag for all three 
examined color-magnitude diagrams (smallest for the (v — 
y)o diagram, and largest for the (b — y)o and /3 diagrams, 
which is expected since the ZAMS in the {v — y)o diagram 
is less steep than in the other two diagrams). 

In Fig. |l] the field stars are plotted together with the 
Pleiades stars adopting the Hipparcos mean parallax of 
the cluster. This plot clearly shows the discrepancy be- 
tween the locus of the field stars and that of the Pleiades. 
The plots in all three color-magnitude diagrams look quite 
similar. In addition to the deviating GO star there are two 
Pleiades stars (Hz II 1338 and 1912) which are ~ 0.6 mag 
brighter than the mean relatio n def ined by the others. 
According to Mermilliod et al. ( 1992|) , one of the stars is 



D. Stello and P. E. Nissen: The problem of the Pleiades distance 



5 




Fig. 1. Color-magnitude diagram of the selected field stars with —0.10 < [Fe/H] < 0.12 and CT(7r)/7r < 0.05 (empty 
circles), compared to the Pleiades (filled stars) adopting the Hipparcos inferred distance modulus of 5.36 mag (from 
Robichon et al. |1999| ). The Pleiades star at the cooler end of the diagram is a GO star and most likely a non-member, 
and the two stars (Hz II 1338 and 1912) about 0.6 mag brighter than the mean relation defined by the others Pleiades 
stars are binaries. 



a spectroscopic binary, and the other a visual binary, in 
agreement with their shift in brightness. 

2.3.1. The stellar magnitude distribution 

To get an estimate of the distance modulus of the Pleiades 
we have to analyze the magnitude distribution, relative to 
the Pleiades MS, of the field stars which define our zero 
point of the ZAMS locus. The observed magnitude dis- 
tribution is a convolution of the underlying evolutionary 
distribution of the field star sample and a Gaussian profile 
coming from the observational errors in My^ . 

A robust least squares fit (second-order) (Freudenreich 

1999| ) of the Pleiades MS is made in the Vq vs. color dia- 



grams. The derived polynomial of the Pleiades MS is used 
to define the shape of the ZAMSpioid- To find the distance 
modulus of the Pleiades the polynomial is shifted by the 
magnitude which make it fit to the ZAMSpicid- The ad- 
vantage of this method is that the robust least squares fit 
of the Pleiades MS will not be significantly affected by the 
stars which lie far from the Pleiades MS. This means that 
e.g. double stars will not bias the locus of the Pleiades MS 
significantly. The disadvantage is that the shifted Pleiades 
MS fit may not match the ZAMSpicid perfectly, because 
the shape of the fitted Pleiades MS is sensitive to the small 
sample of Pleiades data points (29 stars). The shifting of 
the Pleiades MS is done by calculating the individual dis- 
tance moduli of every field star with respect to the poly- 
nomial fit of the Pleiades MS. This give the magnitude 



distribution of the field stars relative to the Pleiades MS. 
All the distance moduli are evaluated as input in the likeli- 
hood function which is the simultaneous probability func- 
tion of all data points. The probabilities of the individual 
data points are described by the result of the convolution 
mentioned above. The underlying evolutionary distribu- 
tion is approximated by oc exp^—x/r) with a sharp edge 
at the ZAMSpioid locus, where t is the fall-off rate due to 
evolution and binarity of the star sample. The width of 
the Gaussian is characterized by the observational error 
a{Mva )- The mathematical expression of the convoluted 
function is: 



f{u) = — exp(--3 ) 

^^2 



Pi 
P2 
Pi 

2p, 



pI 



U-P3, 



1 



■ expi 



A 

2pl 



P2 
U-P3 
P2 



2np4 



cxp( 7T)dv 

^ 2pr 



)(1 + Errf( 



V2p, 



■)) (8) 



where y = u — p3 — p|/p2, Pi is the normalization con- 
stant, P2 = T , p3 is the ZAMSpicid locus, p4 = a{MvQ), 
and Errf(t) is the IDL error function. A best fit is ob- 
tained when the likelihood function takes its maximum 
value which we find by changing r, a{Mvo), and the locus 
of the sharp edge of the ZAMS Field as free parameters. 

The underlying evolutionary distribution (dashed line) 
together with the convoluted function of best fit (thick 
line) is shown in Fig. |^. Additionally, the individual dis- 
tance moduli are binned and over plotted to illustrate the 
field star magnitude distribution for comparison (thin line 
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Fig. 2. The field star magnitude distribution, with re- 
spect to the polynomial fit of the Pleiades MS, in the Mvq 
vs. (& — y)o diagram. Empty circles shows the individual 
distance moduli binned, with a bin size of 0.1 mag and 
connected with the thin line. Dashed line indicates the as- 
sumed underlying star distribution without observational 
errors defining the point of the ZAMSpioid- The thick line 
is the best fit to the data of the convolution of the un- 
derlying star distribution with a Gaussian profile which 
corresponds to the observational error. 

and empty circles). The point of the ZAMSpidd, is the ver- 
tical dashed line. 

The fitted (T{Mvo)Fic\d for all three color-magnitude 
diagrams end up in the range 0.12-0.13 mag which is 
near the expected errors derived in Sect. 2.3 on page 



^. It indicates that the shape of the Pleiades polyno- 
mial match the ZAMSpicid shape very well in all three 
color-magnitude diagrams. The Pleiades distance modu- 
lus found by the fits are: 5.62 ± 0.02 mag, 5.61 ± 0.02 
mag, and 5.61 ± 0.03 mag in the {b — j/)o, {v — y)o, 
and (3 diagram respectively. The stated errors are the 
quadrature sum of two errors. The first coming from 
the uncertainty in the vertical positions of the Pleiades 
MS fits, which is the star to star scatter around the 
Pleiades MS fit divided by the square root of the num- 
ber of stars (5'vo,pioiadcs(Pleiades MS fit)/^/#stars). For 
the three color-magnitude diagrams these errors are: 0.019 
mag, 0.017 mag, and 0.031 mag in the (6 — j/)o, {v — y)o, 
and j3 diagram respectively^. The second error is found 
as the points where the log likelihood function has fallen 
by 0.5 from its maximum value by changing the Pleiades 
distance modulus step by step around its optimum value, 
and optimizing the two other parameters for every step. 
The ±1(7 interval is approximated by a symmetric interval 
around the maximum value by a parabola fit of the log 



^ These errors are from slightly around one to 3/2 times 
larger than the pure observational scatter in Vo,pieiades, and 
show different real effects like binarity, which introduces extra 
scatter around the perfect ZAMS locus. 



likelihood function. This approximation is less than the 
resolution. For the three color-magnitude diagrams these 
errors are: 0.014 mag, 0.015 mag, and 0.015 mag in the 
{b — y)oi (v — y)o, and /3 diagram respectively. 

It is important that the left tail of the field star's mag- 
nitude distribution is fitted well because this is the region 
where the ZAMSpicid has to be found. The extremely good 
fit at the left tail of Fig. |^ supports the trustfulness of the 
method used in this investigation. The right tail is more 
affected by the evolutionary and binary distribution of 
the field star sample, and it could be argued that the as- 
sumed exponential fall-off at the right part of the profile 
is rather simplified, but changing it would not affect the 
fitted ZAMSpicid locus significantly because this part of 
the diagram is separated by several standard errors from 
the ZAMSpioid locus. 

The method of finding the distance modulus presented 
here seems very robust because the distance moduli, errors 
in the distance moduli and a{Mvo) are consistent with one 
another and their estimates in the three color-magnitude 
diagrams. By changing the underlying evolutionary distri- 
bution it is seen that the estimated distance moduli are 
quite stable. We consider the following different underly- 
ing star distributions. Assume the underlying fall-off to 
be in two steps. First a rather steep fall-off followed by a 
less steep fall-off. This scenario would fit the Gaussian to 
be wider (but could still be consistent with the estimates 
given in Sect. 2.2 on page and the Pleiades distance 
modulus would be fitted to be slightly larger (ZAMSpieid 
closer to observed maximum density). The only way to 
get a lower estimated distance modulus is by assuming a 
slower fall-off for the right tail of the underlying distribu- 
tion; the extreme being a wide box function with one edge 
at the ZAMSpioid- Under this assumption the ZAMSpieid 
will be at the point of half maximum of the observed distri- 
bution, which in this case means a lower distance modulus 
by less than 0.05 mag, but at the cost of an unaccept- 
able bad fit to the data. One can of course get an even 
lower estimate if it is assumed that the underlying star 
density will increase on the right side of the ZAMSpieid 
in Fig. 1^, which means that the maximum density of the 
underlying distribution is above the ZAMSpidd- But the 
extremely good fits at the left tail (Fig. |2|) and the con- 
sistent determinations of the Gaussian widths tell us that 
the underlying star distribution must have a sharp edge 
as indicated in Figure |2| To see if the field star sample 
does indeed contain stars not evolved significantly away 
from the ZAMS, the sample is compared with a series of 
isochrones in Sect. 2.3.2| , and the assumed underlying evo- 
lutionary distribution is tested by the aid of evolutionary 
tracks. 

Figure || shows a plot similar to Fig. ^ but instead 
of the adopted Pleiades distance modulus of Hipparcos 
the estimated value of 5.62 mag, which make the fitted 
Pleiades MS polynomial match the ZAMSpidd, has been 
used. Furthermore, the polynomial fit to the Pleiades MS 
is shown, shifted by 5.62 mag (solid line) and by 5.36 mag 
(dashed line). It is easy to see that the Pleiades shifted by 
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Fig. 3. Color-magnitude diagram of the selected field stars with —0.10 < [Fe/H] < 0.12 and a(7r)/7r < 0.05 (empty 
circles), compared to the Pleiades (filled stars) (excluding the likely non- member) adopting the distance modulus of 
5.62 mag, which make the Pleiades MS match the ZAMSpidd- The solid line is the polynomial fit to the Pleiades MS, 
and the dashed line indicates the location of the Pleiades MS according to the Hipparcos distance. 



5.62 mag gives a much better fit to the field stars compared 
to the Hipparcos value. 

2.3.2. The age range 



The age of the field star sample which is used to deter- 
mine the Pleiades distance modulus is analyzed by plot- 
ting isochrones of different ages together with the field 
stars. 

The isochrones used are taken from Lejeune & Schaerer 
(2001), those they denote as 



basic grid" with solar 
metallicity. The effective temperature coming from the 
isochrones is transformed to th e (b — y)o color index by 
the calibration of Alonso et al. ( 1996 ) using a mean value 
of the ci index in their Eq. (9). The possible systematic 
error in (6 — y)o from this transformation is of the order 
of 0.02 mag. 

Figure ^ shows the field stars in the color-magnitude 
diagram and five isochrones of ages: 100 Myr, 500 Myr, 
1 Gyr, 1.5 Gyr, and 2 Gyr. It is evident that the sam- 
ple of field stars contains many stars with ages below or 
around 1.5 Gyr, which are all the stars at the left side 
of or around the 1.5 Gyr isochrone. Since we safely can 
assume that the age of the field stars is distributed uni- 
formly in the plotted color range there must also be many 
field stars of ages around or less than 1.5 Gyr at the right 
side of the plot say in the range 0.30 < (6 — y)o < 0.40. 
There may be a systematic offset of the isochrones with 
respect to the ZAMSficU due to systematic errors in the 
Teff calibration, but Fig. ^ shows that the evolutionary ef- 




0.25 



0.30 
(b-y)o 



0.35 



0.40 



Fig. 4. Color-magnitude diagram of the selected field 
stars with -0.10 < [Fe/H] < 0.12 and cr(7r)/7r < 0.05, 
compared to five isochrones representing the ages: 100 
Myr, 500 Myr, 1 Gyr, 1.5 Gyr, and 2 Gyr. 



fects on the isochrones from 100 Myr (approximate age of 
the Pleiades cluster) to 1.5 Gyr is rather negligible in the 
color range 0.30 < {b — y)o < 0.40. In that range we still 
see the discrepancy between the position of the Pleiades 
and the ZAMSpieid if the Hipparcos distance is adopted. 
Hence, we conclude that the discrepancy cannot be ex- 
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plained as due to evolution away from the ZAMS of the 
field star population. 

As an additional check, evolutionary tracks from 



Lejeune & Schaerer ( |200l| ) have been used to calculate 
the theoretical stellar magnitude distribution at a given 
{b — y)o in order to see how well it agrees with the un- 
derlying evolutionary distribution of My adopted in Sect. 
2.3.l[ We assume a constant star formation rate for so- 
lar metallicity stars over the lifetime of the galactic disk 
(~ 8 Gyr) in agreement with the age-metallicity diagram 
(Fig. 14) of Edvardsson et al. (1993), and a constant initial 
mass function over the small mass range corresponding to 
a given (b — y)o. The calculated magnitude distribution 
is similar to the underlying evolutionary distribution of 



My adopted in Sect. 2.3.1, i.e. with a sharp edge at the 
ZAMS and a steep evolutionary fall-off. The fall-off is less 
steep at the blue end of the {b — y)o range and somewhat 
steeper at the red end, and the fall-off has a tendency 
of a two-step function; first a steepe r part followed by a 
less steep part. As discussed in Sect. 2.3.1 on page ^ this 
could indicate that our fitted Pleiades distance modulus is 
slightly underestimated. Altogether, we conclude that the 
assumed underlying evolutionary magnitude distribution 
of our field stars is supported by models for the stellar 
evolution. 



2 A. [Fe/H] sensitivity of ZAMSpieid 

It is known that the ZAMS locus is dependent on the 
metallicity, such that low metallicity stars define a fainter 
ZAMS than the high metallicity stars. To be able to 
test if the discrepancies in the Pleiades distance could 
be explained as a possible error in the adopted Pleiades 
metallicity, we have analyzed how much the locus of the 
ZAMSpicid changes as a function of the metallicity. 

The selected sample for this part of the investigation 
consis ts of all F-type field stars in the catalogue of Olsen 
( |1999D , which have a relative error in the parallax mea- 
surement less than 5%. This selection gives a sample of 
2309 stars. Five plots, each representing field stars in dif- 
ferent metallicity intervals, are then made. The intervals 
are ±0.10 dex wide in [Fe/H], ranging from —0.45 to -1-0.15 
dex, and with a 0.10 dex overlap from one interval to the 
next. 

Distance moduli for the Pleiades are found for the five 
metallicity intervals in all three color-magnitude diagrams. 



The method is similar to that presented in Sect. 2.3.1 
and the individual distance moduli from every diagram 
is plotted together in Fig. ^. The indicated error bars in 
Fig. ^ are larger for the data points representing the lower 
metallicity intervals, which is a result of less stars falling 
in the lower metallicity bins. By analyzing the fits to the 
field star density profile of the lower metallicity inter vals 
it was observed that the fitting function (see Sec 2.3.1 Eq. 
(^) did not fit the left tail of the distribution as good as 
shown in Fig. § simply because of the few data points avail- 
able. Linear least squares fits are made to the result of all 
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Fig. 5. Metallicity sensitivity of the ZAMSpioid in the 
three investigated color- magnitude diagrams: (3 (dia- 
monds), {b — y) (squares), and {v — y) (X's). Solid and 
dashed lines are least squares fits (see text). 

three color-magnitude diagrams in the metallicity range 
-0.35 < [Fe/H] < 0.05 (dashed hues), and additional fits 
where the two lowest metallicity points are ignored (solid 
lines). 

The metallicity sensitivity from the three color- 
magnitude diagrams is found to be (solid lines): 



^ (^^Vo)zAMS 
^ (^\^o)zAMS 
^ (^Vo)zAMS 



-0.01(9) A[Fe/H], from (3 
-0.40(9) A[Fe/H], (6 - y)o 
-0.76(9) AiFe/H], and {v - y)o. 



The trend is the same for both solid and dashed lines, 
though solid lines show a higher metallicity sensitivity in 
{b ~ y)o and {v — j/)o and no significant sensitivity in the 
(3 diagram. We believe that the relations indicated by the 
solid lines are the most correct since the lower metallicity 
bins may suffer from small evolutionary effects due to a 
larger percentage of more evolved stars. That the highest 
metallicity sensitivity is observed in the (v — y)o diagram 
and the lowest in the /? diagram is expected since line blan- 
keting from the metal lines affects the spectrum more at 
the short wavelengths, thus v is more affected than b which 
is again more affected than y. We note that Grenon ( 1999[ ) 
finds a much larger shift A {Mva)zAMS " 

-1.67A[Fe/H] 

based on Geneva photometry, though the error of the 
metallicity sensitivity coefficient is not given. His relation 
refers, however, to early K-dwarfs and it is not clear on 
which color index it is based. A direct comparison can 
therefore not be made to the shifts we are finding for the 
F-type stars. 

Figure || supports that the metallicity of the Pleiades 
is around the solar value as determined from the mi index 
if the adopted Hyades metallicity is ^ 0.12 dex. Only for 
this metallicity can we get consistent results of the dis- 
tance modulus from all three color-magnitude diagrams. 
The inferred distance modulus is in the range 5.58-5.64 
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mag. We also see that a Hipparcos distance modulus of 
5.36 mag is not consistent with one metallicity value. The 
{v — y)o diagram demands [Fe/H] ~ —0.35, the (b — y)o 
diagram [Fe/H] ~ —0.70, and the /? diagram can not be 
fitted to a distance modulus of 5.36 mag no matter the 
chosen metallicity of the field stars. 

From the above metallicity sensitivity analysis one 
could conclude that the /3 diagram gives the most reliable 
estimate of the distance modulus, because it is rather in- 
sensitive to metallicity deviations between field stars and 
Pleiades stars, and furthermore it is insensitive to inter- 
stellar reddening. It should, however, be remembered that 
the P diagram has the largest error in the distance modu- 
lus determination (0.03 mag); 3/2 the size of the {b — y)o 
and {v — y)o diagram errors. If we combine the distance 
moduli from all three color-magnitude diagrams the final 
distance modulus is: (m — M)o = 5.61 ± 0.03 mag|^ 

3. Discussion and conclusion 

The extensive multi color MS fitting analysis based on 
Stromgren photometry (this paper) gives a Pleiades dis- 
tance modulus of 5.61 ± 0.03 mag (the mean value from 
all three color-magnitude diagrams), which is in very good 
agreement with the distance modulus given by the multi 
color MS fitting analysis of Pinsonneault et al. ( 199^ ) who 
find a distance modulus of 5.60 ± 0.05 mag. On the other 
hand, the distance modulus derived from the Hipparcos 
parallaxes are all in the range of 5.33-5.37 mag, which is 
not consistent with the former results, and the quoted er- 
rors. Even if the quoted errors from MS fitting are not 
representative for the actual uncertainty in this method, 
it must, from the current investigation, be concluded that 
the Hipparcos distance modulus is not consistent with the 
ZAMS of nearby field stars (Fig. |l|). Thus, either are the 
Hipparcos parallaxes affected by systematic errors, which 
are significantly larger than expected, or there is some- 
thing unexpected about the Pleiades cluster. 

If the discrepancy is caused by some anomaly of 
the Pleiades, one possible explanation is the metallicity. 
Extensive investigation of this possibility has been per- 
formed (e.g. Pinsonneault et al. 



be ruled out. Hence, we cannot confirm the recent claim 



1998 and the current in- 



vestigation presented in Sect. |2^), and from these results, 
supported by the spectroscopic metallicity determination 
by Boesgaard & Friel ( |1990D , it seems very unlikely that 
the adopted Pleiades metallicity around the solar value is 
more than 0.05 dex from the true value. Even a deviation 
of 0.1 dex is far too small to explain the discrepancy and it 
would lead to inconsistent results from the different color- 
magnitude diagrams. In addition, the relative comparison 
between the Pleiades and field stars of the same metallicity 
(Sect. ^ ensures that the possible difference in the relative 
metallicities are so small that a metallicity deviation can 



by Grenon (1999) that the Pleiades problem is solved by 
adopting a low metallicity, [Fe/H]pioiadcs = —0.11, of the 
cluster. Furthermore, the large metallicity difference be- 
tween the Hyades and the Pleiades based on Geneva pho- 
tometry, A[Fe/H] ~ 0.25 ± 0.03, is inconsistent with the 
difference, A[Fe/H] = 0.11 ± 0.03, which we derive from 
Stromgren photometry. 

Another possibility is an abnormal helium abundance 
of the Pleiades. To see which value of Y would be required 
to change the Pleiades MS locus by 0. 3 mag, a calculation 
is made by Pinsonneault et al. ( |l998|) , who find th at the 
value is as high as y ~ 0.37. The study of Nissen ( 1974 ) 



revealed no intrinsic scatter in Y greater than approxi- 
mately 10% in nearby MS field B stars; much smaller than 
the 30%-40% change in Y required for the Pleiades. There 
are, however, investigations which indicate large cluster 



to cluster s catter in the helium abundance (Nissen 1976 
Lyubimkov 1977| ), and it has been suggested that this is 



The distance modulus found by CP76, based on all 30 mem- 
ber stars in their Table II, was (m — M)o = 5.53±0.04 mag, and 
by rejecting four stars, probably double stars or non-members, 
their result was (m — M)o ~ 5.60 ± 0.03. 



the key explanation to the Hyades ci-anomaly (Stromgren 
et al. 1982). To test if this is the explanation of the 
Pleiades problem too, an attempt should be made to mea- 
sure the surface helium abundance of the hot stars in the 
Pleiades and other young cl uster s spectroscopically. 

Recently, van Leeuwen (1999) has suggested that the 
Pleiades problem is caused by an age effect, and claims 
that other very young open clusters show the same devi- 
ation as the Pleiades. The investigation of van Leeuwen 
( |1999D is based on a comparison in the color-magnitude 
diagram (with the metallicity-sensitive B ~ V color) of 
nine open clusters, all shifted to their Hipparcos mean dis- 
tance. But this is done without correcting for differences 
in the metallicity abundances fi rst. A test of the age effect 
is made by Pinsonneault et al. ( ^000 ), based on 8 clusters 
(incl. the Hyades and Pleiades) and no age effect is seen 
in the difference between the MS fitting and Hipparcos 
distances. If the youth of the Pleiades has a significant af- 
fect on the Pleiades MS locus, one might also expect that 
young field stars would show the same effect. The study of 
Soderblom et al. (1998) of chromospherically active (and 
therefore assumed young) stars gave, however, no indica- 
tions of that. In our own sample of field stars very few 
are expected to be as young as the Pleiades so we cannot 
test the suggestion of van Leeuwen ( 1999| ). It should be 
emphasized, however, that the explanation given by van 
Leeuwen (1999) of the Pleiades problem as an age effect, 
mostly relies on the relative shift between the Pleiades 
and Hyades in the temperature range where no obvious 
evolutionary effects away from the ZAMS are seen (cor- 
responding to (6 — y)o > 0.30 mag or (J5 — V)q > 0.50 
mag). Furthermore, the age explanation of van Leeuwen 
( p^99| ) is in disagreement with theoretical models of stel- 
lar evolution, which predict negligible evolution away from 
the ZAMS during the first couple of billion years of the 
lifetime of late F-type stars (see isochrones in Fig. ^) . 

Some investigations give hints of spatial systematic er- 
rors in the Hipparcos parallaxes which are larger than 
expected (Pinsonneault et al. 1998; Narayanan & Gould 
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199£). In addition, there is a statistical correlation caused 
by the imperfect distribution of data points over the el- 
lipse described by the parallactic motion. In particular, 
for a star on the ecliptic, equal numbers of measurements 
should be obtained on both sides of the Sun. This was 
not fulfilled for Hipparcos, and caused correlations (p^) 
between right ascension and parallax (ESA 1997, Vol.1, 
p. 325). In Pinsonneault et al. (1998) it is questioned if 
this type of correlation would have an effect on the paral- 
lax values. But from the extensive test by Robichon et 
al. (1999) of this issue it must be concluded that the 



correlations do not introduce significant errors in the 
Hipparcos parallaxes. 

What if the whole controversy, about the Pleiades dis- 
tance modulus, is caused by a real effect? What will the 
effects on the MS fitting and Hipparcos results be, if the 
cluster is non-spherical (sphericity has until now been im- 
plicitly assumed)? The angular size of the Pleiades clus- 
ter can be approximated by the size of the region cov- 
ered by the member stars used in the investigation of e.g. 
Narayanan & Gould ( 1999 ). This gives a radius of ^ 6° 
which, at a distance of 130 pc, corresponds to a radius of 
~ 14 pc. A typical real difference in the cluster member 
distances would then be around 14 pc which corresponds 
to a 1 mas difference in the parallax. There is no rea- 
son why non-sphericity of the Pleiades should not be the 
case. There is plenty of evidence that open clusters can be 
non-spherical e.g NGC 2264 (the Cone Nebulae), and ac- 
tually Raboud & Mermilliod ( 1998D have shown that the 
distribution of the Pleiades stars projected on the sky is 
elliptical with an ellipticity of 0.17. It could be that the 
cluster has a more oblong shape in the direction of the line 
of sight say with a length that is twice the projected diam- 
eter. One could then imagine the following scenario: The 
first born bright stars (O and B-typc) forms in one part 
of the gas cloud, and they start to blow the gas cloud in 
one initial direction, and therefore these stars will end up 
at one end of this deform shape (as observed in the NGC 
2264 case) , and the fainter stars (F and G-type) will form 
a "tail" (as an overall trend). So if we see this shape head- 
on there will be a trend that the brightest B-type stars 
are closer to us, and the later classes are further away. 
Because the calculation of the Hipparcos mean parallax 
gives the largest weight to the brighter stars, the result 
will be a slightly shorter distance, than the actual mean 
cluster distance (Pinsonneault et al. 1998, Fig. 20). On 
the other hand, the MS fitting method relies mostly on 
the fainter stars (A to G-type), which are farther away. 
These stars are located in the color-magnitude diagram, 
where the slope of the ZAMS is less steep, and therefore 
they give the smallest errors in the distance estimates 
(additionally for many clusters the hotter stars are also 
evolved away from the ZAMS). So this means that the 
distance found from MS fitting will be larger than found 
from Hipparcos parallaxes, and possibly closer to the ac- 
tual mean cluster distance. Such oblong shape of the clus- 
ter, could indicate a kinematic history that does not follow 
the usual assumptions, which also explains why Robichon 



et al. ( 1999| ) find unusual features in the kinematics of the 
Pleiades. The consequence of this is that the assumptions 
used by Narayanan & Gould ( 1999D docs not hold. One 
could argue that if the Pleiades have such non-spherical 
shape, it would be expected that some fainter stars (which 
have larger distances) were present near the center of the 
cluster in e.g. Fig. 20 of Pinsonneault et al. (199J ). But if 
there is a dark cloud just behind the bright stars as in the 
NGC 2264 case, none of these stars will be observed. 

The idea of a deform Pleiades cluster is a tempting 
answer to the Pleiades problem, because it includes most 
of the evidence presented in this discussion. What is de- 
scribed above, as a possible non-spherical Pleiades clus- 
ter, shall be viewed as an illustration or example of the 
possible effects on the distance determination, due to a 
deform and non-symmetrical cluster. Though there are 
studies of the velocity dispersion among Pleiades clus- 
ter members (van Leeuwen 1994) and mass segregation 
(Raboud & Mermilliod 1998) suggesting that the Pleiades 
is a bound and quite relaxed system, it could be very inter- 
esting to investigate the possibility of a non-symmetrical 
cluster by e.g. extensive kinematic analysis of the Pleiades. 
Furthermore, future astrometric space programs will be 
capable of determining the distances to the individual 
Pleiades stars with an improved accuracy of 2 to 3 orders 
of magnitude, compared with Hipparcos. These measure- 
ments will provide a very good three-dimensional picture 
of the Pleiades cluster. 
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